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We discuss the construction of optimized electronic filters 
using inverse scattering methods. We study a wide range 
of densities and temperatures, room temperature included. 
Discretization methods of the potential (including the self- 
consistent potential of the conduction electrons) are worked 
out that retain all its properties. 
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I. INTRODUCTION 

A semiconductor heterostructure can be modeled by 
a system of equations describing (with a certain degree 
of completeness and precision) the state of the system. 
The equations depend on a set of structural and com- 
positional data: the spatial dependence of the chemi- 
cal composition (including dopant profiles), the applied 
external fields, etc. The system behavior (response) is 
described by functional data, such as the electric or ther- 
mal conductance, the energy dependence of the electron 
transmittance, the wavelength dependence of the opti- 
cal absorption coefficient, etc. The functional data can 
be computed using the solution of the equations, and are 
thus functionals of the structural and compositional data. 

To design a heterostructure for a certain application 
is to find a set of structural and compositional data, 
which is physically (and technologically) achievable, such 
that the values of a chosen subset of functional data will 
be within a desirable range. The designer, in principle, 
solves an inverse problem: inverting the dependence of 
the functional data on the structural and compositional 
data. 

In the usual approach, the starting point is a proposed 
configuration of the structural and compositional data 
which is then varied to optimize the component's perfor- 
mance; at each iteration the spectral-scattering data de- 
scribing the performance have to be re-calculated. In this 
method, much depends on the initial choice of parame- 
ters as in general several non equivalent local extrema of 
the response are present in the parameter space. 

Inverse methods, when available, are not subject to 
this drawback. Contrary to the previous situation, the 
designer starts from data describing the system's perfor- 
mance, such as current spectral response curves and de- 
termines the optimal structural and compositional data. 
He automatically gets an optimized set of parameters. 



Even if the model used to make the optimization cal- 
culation feasible lacks all the details of the full physical 
situation, one has still the possibility to use the tradi- 
tional approach to fine tune the results from such a new 
starting point. 

The inverse method we use requires three different 
steps [1,2]: 

1. The reconstruction of the phase of the transmittiv- 
ity of a quantum device using Finite Dispersion Relations 
under the very convenient form of the Pade approxima- 
tions. 

2. The Standard Inverse Scattering method, making 
use of the results of Kay and Sabatier [3,4] for the in- 
verse problem with rational coefficients (Pade approxi- 
mations). This allows to build an efficient algorithm for 
solving the inverse problem and obtaining the potential. 

3. A unitary transformation that maps the place de- 
pendent mass problem of a BenDaniel equation onto the 
usual constant mass Schroedinger problem keeping the 
potential in the equation local [2] . 

By applying the previous procedure, it becomes possi- 
ble to design new realistic quantum wells [5,6] submitted 
to optimization constraints. In particular, one can use 
inverse methods to design and fabricate novel color sensi- 
tive QWIPs in the medium- and long-wave IR. Moreover, 
one of the most unexpected and astonishing results of the 
inverse method is that inclusion of the self-consistent po- 
tential of the conduction electrons is far much simpler 
than in the usual direct approach [2]. 

As a first application, we will use the versa- 
tility of the inverse scattering methods for the 
Schroedinger/BenDaniel-Poisson model [1,2], for design- 
ing and building prototypes of improved electronic filters. 

An aspect that has not yet been dealt with in depth 
is how to discretize the potential profiles obtained by in- 
verse scattering methods so as to make their production 
viable while mantaining at the same time their proper- 
ties, their electron wave reflectance in particular. An 
added complication is the necessary presence in the de- 
vices of substantial densities of conduction electrons: the 
electric potentials they generate do not, if suitably com- 
pensated by doping the devices [2], significantly alter 
the reflectance of the "smooth" potentials obtained by 
inverse scattering; but what about the discretized po- 
tentials? The aim of the present paper is to deal with 
the problems stemming from this necessity to both dis- 
cretize the potentials and have substantial densities of 
conduction electrons. To avoid unnecessary complica- 
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tions, we shall skip the last step of the procedure de- 
scribed in Ref. [2], namely the mapping of the (constant 
mass) Schroedinger equation into the (variable mass) 
BenDaniel and Duke equation, as it has no direct bearing 
on the present discussion. 

Throughout the present paper we shall make use of 
the "natural" units thus defined: energy (c.u.): 0.1 eV; 
mass: 0.1 m e ; length: 1 nm: time: 2.3845 fs; in these 
units % = 2.7604. 



II. INVERSE SCATTERING PROBLEM 

The procedure to build a potential having a given re- 
flectance profile for elecron waves and no bound state 
is described in detail in Ref. [2]; here we wish to stress 
that while standard type II Pade [7] fitting routines can 
be a useful first step in building a Pade approximant, it 
is often the case that the parameters thus found must 
be adjusted, as care must be taken to avoid obtaining a 
potential either impossible to build because exceeding in 
height the available potential difference (for Al c Gai- c As, 
where c is the Aluminum concentration, this value is 
about 3e.w.) or impractical to grow because having a sup- 
port too wide and/or too many peaks. The reconstruc- 
tion of the potential makes use of the poles and relative 
residues of the reflection coefficient R + (k) for plane waves 

incoming from +oo where k = \J mE/% 2 is the electron 
momentum. The frequencies of oscillation relate to the 
real part (aj) of the poles (where j = l,n, n being the 
number of poles), the decay factors to the imaginary part 
(bj), and the amplitudes to the residues (pj). Large val- 
ues of the Oj 's give fast oscillations of the potential; small 
values of the bj's result in slow decays of the oscillations 
and therefore in a wide support for the potential; finally, 
large pj's -that are often the case when poles are very 
close to each other- result in large oscillations of the po- 
tential; all three must therefore be avoided. Sometimes 
these conditions can be satisfied with no problems, some- 
times instead they force us to some kind of compromise. 
An example of the latter case is a bandpass filter: since 
the more the poles are close to the real axis and close 
to each other the more the transmission plateau is flat 
and close to unity, some extra care must be taken when 
choosing the Pade approximant parameters. 



that of the original potential gives us an (almost) zero 
transmittance in the forbidden zones and an almost sat- 
isfactory plateau in Fig. 2 without significant shifts. A 
desirable characteristic of this method is the unexpected 
reduced width of the two peaks in Fig. 1. A second scheme 
can be used to enhance this characteristic: 

b) If we move the position of the steps to the zeroes 
of the original potential, both peaks are now shifted, but 
we have succeeded to reduce the width of both of them by 
a factor of 4! 

As a rule, both methods should be tested and the one 
that gives the result nearest to the desired one retained. 

IV. DOPING FOR NON-ZERO DENSITY OF 
CONDUCTION ELECTRONS 

The (selfconsistent) electric field generated by the con- 
duction electrons can seriously modify the effective po- 
tential and thus alter the transmission profile. It must 
therefore be compensated. We achieve this in two steps: 
first, from the potential obtained by inverse scattering 
we subtract the selfconsistent field of the electrons them- 
selves and of some suitable doping charge such that the 
total charge and dipolc moment in the device are zero 
(see Ref. [2]); then, after obtaining the "chemical com- 
position potential" by discretization of the resulting po- 
tential, we calculate the selfconsistent field for this new 
potential itself, adjusting the position and value of the 
doping charge, so that again the total charge and dipolc 
in the device are zero and the effective potential is as close 
as possible to the " chemical composition potential" itself. 
An iterative selfconsistent calculation starting from the 
doping charge and dipolc moment obtained for the poten- 
tial before discretization converges in the worst case we 
considered in less than ten steps to a precision of five dig- 
its. Note that, once the correct doping has been added, 
the changes to the potential are usually very small; on the 
other hand a wrong doping is bound to seriously deteri- 
orate the transmittance profile. In both of the following 
examples we assumed the doping to be in the form of a 
sheet charge; this approximation is valid as long as the 
doping can be kept in a region more narrow than the 
steps of the potential; should this not be the case, a suit- 
able charge width should be considered and taken into 
account, requiring only minor changes to the codes. 



III. DISCRETIZATION OF THE POTENTIAL 



V. TWO EXAMPLES 



A continuous potential as the one obtained by inverse 
scattering [2] is impractical to build; we must therefore 
discretize it by a sequence as small as possible of flat 
steps. Several options are possible, the choice depending 
on the characteristics of the heterostructure we wish to 
build: 

a) placing the steps at the inflection points of the orig- 
inal potential and choosing their height to be equal to 



To illustrate the above procedure, we now give two 
examples. The starting point for both is a Pade approx- 
imant to the reflectance of the form: 



K{E) 



A 2N (E) 



(A2n(E) + EK,2(N+2){E)) 



^2(N+2)(E) = Uj=l,N+2 (E - Fj) 2 + Sj 



(1) 

(2) 
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A 2N (E) = g • [n i=liAr (E - Ei)] 2 Ei, F h 8j e R+ (3) 

Eq. (1) guarantees that 7^.(0) = 1; it was already given 
in Ref. [2] but there higher degrees for IC 2 (n+2) (E) were 
allowed. We choose to use the lowest possible degree 
ensuring the absence of (5-function terms in the potential 
[8]. Eq. (2) guarantees that, excepted the zero at the 
origin, the transmittance has complex conjugate zeroes; 
Eq. (3) instead guarantees that the reflectance has real 
double zeroes. The three equations together guarantee 
that < K(E) < 1 for E > 0. 

We assume the devices to be grown from Al c Ga\^ c As; 
we therefore use as parameters those for the substrate 
GaAs: effective electron mass = 0.67 in natural 
units, and electrical permeability e = 13. 

A. A Filter with Two Narrow Transmission 
Resonances 

Our construction of a filter with two narrow peaks, 
the one at higher energy being the narrower [9], starts 
from a Pade approximant having the same number of 
poles and zeroes as the one given in Refs. [2,1] but in 
the form (1) given above with: N = 2, g = 1.5, E\ — 
OAe.u., E 2 = l.le.u., F 1 = 0.2e.u., <5i = 0.05e.w., F 2 = 
0.8e.u., S 2 = O.Ole.u., F 3 = 1.2e.u., S 3 = O.lOe.u., F 4 = 
1.5e.u., and £4 = 0.20e.u.. E\ and E 2 give the positions 
of the transmittance peaks; the four Fj frame them and 
determine the peaks' width. 

The final transmittance using procedure b) to dis- 
cretize the potential, is shown in Fig. 1 for three tem- 
peratures (T = 70, 233, and 300 K) and two background 
electron densities (10 13 , 10 15 ); the six curves are undis- 
tinguishable. The compensating doping charges and their 
positions are given in table 1. We have also studied 
other temperatures and densities that are not reported 
for lack of space. Only at the highest density studied 
(10 17 el /cm 3 ) the reflectance is strongly affected by the 
electron density, especially at low temperatures. The rea- 
son is that since at low temperatures the deviation from 
the background value of the electron density is less spread 
out, it has sharper peaks which are reflected in analo- 
gous sharp peaks in the selfconsistent electrical field; the 
plateaus of the discretized potential are therefore more 
strongly deformed than at higher temperatures. 

B. A Bandpass Filter 

For this case some preliminary remarks are in order: 
the reflectance of an ideal bandpass filter should be zero 
inside the band and one everywhere else. On the other 
hand for a (piecewise) continuous potential that decays 
at infinity we have: 

a) 11(E) -> for E -> 00. 



As a consequence, the forbidden region will have to be 
finite; it will be our aim to make it as wide as possible 
on both sides of the transmission band. 

b) 11(E) = 1 only for E = 0; otherwise 11(E) < 1. 
Therefore, there will be a residual transmittance in the 

forbidden region; it will be our aim to make it as small 
as possible. 

c) 1Z(E) can be zero only on a set of points which is 
at most countable. 

Because of this, the desired transmission band will 
have to be built by a superposition of peaks, close enough 
that the transmission plateau is almost flat and sharp 
enough that the transition from forbidden region to al- 
lowed band is narrow compared to the band width. 

We found convenient to choose the following param- 
eters in eq. (1): N = 3, g = 20, E 1 = 0.9e.u., E 2 = 
l.le.u., E 3 = 1.3e.u., F 1 = F 2 = 0.3e.it., <5i = 0.08e.it., 
S 2 = 0.22e.u., F 3 = F 4 = F 5 = 1.6e.u., S 3 = 0.48e.it., 
S4 = 0.30e.u., and 65 = 0.16e.u.. The transmittance 
plateau is given by the three reflectance zeroes in Ei, 
i = 1,2, 3; the five Fj frame them and determine the two 
forbidden zones at each side of the allowed band. 

The final transmittance using procedure a) to dis- 
cretize the potential, shown in Fig. 2 for the same 
three temperatures and two background electron densi- 
ties used in the previous example, is again equal in all 
cases; the compensating doping charges and their posi- 
tions are given in table 2. Again we found that only at 
the highest density studied (10 17 el /cm 3 ) the reflectance 
is strongly affected by the electron density. 



VI. CONCLUSIONS AND 
AKNOWLEDGEMENTS 

We have shown that both discretization of the poten- 
tial and substantial densities of conduction electrons can 
be dealt with in ways that preserve the essential charac- 
ters of the desired reflectance. We should note that in 
both the examples given the background electron distri- 
bution is repelled from the potential and therefore the 
required compensating doping is of the N-typc; this can 
sometimes be impractical. A way to circumvent the prob- 
lem is to add a (monodimensional) bound state via a 
Darboux transformation [8]: the high electron density in 
it ensures that the doping be of the P-type. 

Finally, we would be glad to receive feed-back from 
readers. We are ready to analyze and design electronic 
filters that are difficult or impossible to deal with by the 
ordinary method. 

We wish to aknowledge the participation to the first 
phases of the present research of the late G.A. Mez- 
incescu, who has been for a long time a driving force 
in our group. 
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T(K) el. density (el /cm 3 ) charge(cm 2 ) position(nm) 



70 


1 nl3 

10 


- 3.1729 * 10 


12.8 


70 


10 15 


- 3.1455* 10 9 


12.9 


233 


10 13 


- 1.8134* 10 7 


10.4 


233 


10 15 


- 1.8113* 10 9 


10.3 


300 


10 13 


- 1.6038* 10 7 


10.1 


300 


10 15 


- 1.6061 * 10 9 


10.1 



TABLE II. Doping charge and position at different tem- 
peratures and background electron densities for the bandpass 
filter described in the text. 



FIG. 1. The reflectance of our double peak filter at 
T = 70, 233, and 300 K for two different densities of con- 
duction electrons (10 13 el /cm 3 and at 10 15 el /cm 3 ): all six 
curves are indistinguishable on this scale from the "optimal" 
reflectance obtained using procedure a) described in the text 
and truncating the support at 50 nm. The width ratio of the 
two peaks is approximately 3. 



FIG. 2. The reflectance of our bandpass filter at 
T — 70, 233, and 300 K for two different densities of con- 
duction electrons (10 13 el /cm 3 and at 10 15 el /cm 3 ): all six 
curves are indistinguishable on this scale from the "optimal" 
reflectance obtained using procedure b) described in the text 
and truncating the support at 100 nm. 



T(K) el.density(el/cm 3 ) charge(cm 2 ) position(nm) 



70 


10 13 


- 2.3239* 10 7 


9.0 


70 


10 15 


- 2.3008* 10 9 


9.0 


233 


10 13 


- 1.8134* 10 7 


6.8 


233 


10 15 


- 1.8113* 10 9 


6.8 


300 


10 13 


- 1.6038* 10 7 


6.5 


300 


10 15 


- 1.6061 * 10 9 


6.6 



TABLE I. Doping charge and position at different temper- 
atures and background electron densities for the filter with 
two narrow transmission peaks described in the text. 
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energy (100 meV) 




energy (100 meV) 



